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ABSTRACT 



We combined sensitive near-infrared data obtained with ground-based imagers 
on the ESO NTT and VLT telescopes with space mid-infrared data acquired with 
the IRAC imager on the Spitzer Space Telescope to calculate the extinction law 
A\/A Ks as a function of A between 1.25 and 7.76 /xm to an unprecedented depth 
in Barnard 59, a star forming, dense core located in the Pipe Nebula. The 
ratios A\/Ak s were calculated from the slopes of the distributions of sources 
in color-color diagrams A — Kg vs. H — K$- The distributions in the color- 
color diagrams are fit well with single slopes to extinction levels of Ak s ~ 7 
{Ay pa 59 mag). Consequently, there appears to be no significant variation of the 
extinction law with depth through the B59 line of sight. However, when slopes are 
translated into the relative extinction coefficients A\/Ak s , we find an extinction 
law which departs from the simple extrapolation of the near-infrared power law 
extinction curve, and agrees more closely with a dust extinction model for a 
cloud with a total to selective absorption i?y=5.5 and a grain size distribution 
favoring larger grains than those in the diffuse ISM. Thus, the difference we 
observe could be possibly due to the effect of grain growth in denser regions. 



Finally, t 
study of 



re slopes in our d i agram s are somewhat less steep than those from the 



Indebetouw et al 



(120051 ) for clouds with lower column densities, and 
this indicates that the extinction law between 3 and 8 /zm might vary slightly as 
a function of environment. 



Subject headings: dust, extinction — infrared: ISM 
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Introduction 



One common technique to study the structure of molecular clouds is to map their 
distribution of dust extinction (which in turn is proportional to the column density of 
H 2 ) by measuring the reddening of background stars. This reddening is observable as 
an excess in a given photometric color, which is translated to extinct ion by assuming 



a conversion factor from a previously determined extinction law (e.g. 



Lada et al. 



1994 



Lombardi fc Alves 



20011 ). Towards many molecular clouds dust extinction is so large that 



reddened background stars are only detectable in the infrared, making it necessary to image 
clouds with specific detectors for that regime. Under the best conditions (high data quality, 
high density of background sources based on the location of the field against the Galactic 
plane, etc.), it is possible to use near- infrared (NIR)S cameras on large aperture telescopes 



(8m class) to observe colors of stars in the b ackground of mo 



extinctions reaching up to A v =50 mag (e.g. 



Alves et al. 



ecular cores with maximum 



20021 ) . Higher extinction regimes 



(50< Ay <100 mag), such as those found in so me very dense, proto stellar cores and more 



Teyssier et al. 



20021 ). can not be probed 



recently in Infrared Dark Clouds (IRDCs) (e.g. 
even with those instruments and the best observational conditions. If one wishes to measure 
extinction at those levels it is necessary to consider the use of MIR wavelengths, which are 
capable of penetrating deeper into the cores. 

The Infrared Array Camera (IRAC) on board the Spitzer Space Telescope, with its 
relatively large field of view (5.1 x 5.1') and its capability of simultaneous observation of 
four MIR bands ([3.6], [4.5], [5.8] and [8.0] /im) allows us to detect those deeply extincted 
sources that are otherwise hidden in the optical and NIR, and to measure with reasonable 



Tn this paper we define Near-infrared (NIR) as the range from 1.2 to 2.2 /zm (J,H,K) typ- 
ically covered from ground, and Mid-infrared (MIR) as the range from 3.5 to 8.0 /zm covered 
by IRAC. 
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precision color excesses that correspond to the highest extinction levels of the most opaque 
pre-stellar and young star forming cores. Combined NIR+MIR colors of the form A — K 
can be used to determine the color excesses of those sources with no J (1.2 /im) or H 
(1.6 /im) detection, while stand alone MIR colors can be used for sources hidden even at 
K (2.2 /im). The conversion of color excesses to extinction for deeply embedded sources, 
however, requires of knowledge of the extinction law in the NIR and MIR. 



e.g. 



n the NIR regime, the extinction law has been determined in nu merous studies 



Rieke Lebofsky 



1985 



Kenyon et al. 



1998 : 



Lombardi et al. 



20061 ). and it has been 



shown to be an approximately uniform power law of the type with j3 ~ —1.8. The 
NIR extinction law appears to be universal and independent of Ry, the total to selective 



extinction ratio, although slight variation s from place to place and as a function o 



(i.e. extinction level) have been reported (IMoore et al. 



2005 



Nishiyama et al. 



2006|)- 



depth 



In the MIR regime, the determination of the extinction law has been more difficult 
because this wavelength range is best accessed from space. Spectroscopic studies with 
the SWS spectrometer on the Infrared Space Observatory (ISO) have been used to 
determine the extinction law by comparing observed and expected intensities at specific 



MIR w avelengths for different dust extinction models. In one of these studies, iLutz et al. 



(1996 



, hereafter L96) obtained the 2.5-45 /im spectrum of SgrA*, and their observed 
extinction law was found to be almost flat between 4 and 8 /im, instead of continuing 
the NIR power law towards longer wavelengths, as expected from models of clouds with 
total to selective extinction Ry=3.1 (typical o f the diffuse ISM) and a standard mixture of 



graphite-silicate grains (e.g. 



Li fc Draine 



20011 ) . The extinction law of L96 instead displayed 



a better agreement with models which are constructed with both a smaller number of very 
small grains and a more significan t contribution of large grain s as might be expected for 



more dense clouds (Ry > 4.0) (see 



Weingartner &: Draine 



2001 



, hereafter WD01). However, 



the spectroscopic studies by iBertoldi et al.l ([1999) and iRosenthal et al.l (120001 ) produced 



significantly different results: they calculated relative line intensities of H 2 emission towards 
the Orion Molecular Cloud (OMC) also from SWS spectra, but their measurements are in 
better agreement with an extinction law that extends the power-law observed in the NIR, 
consistent with diffuse medium models. 



Indebetouw et al.l (120051 . hereafter 105) combined IRAC and 2MASS for two distinct 
fields of the Galactic plane to determine the MIR extinction law photometrically from color 
excesses of background stars as done in the NIR. They derived the extinction law in the 
form A\/Ak as a function of A between 1.2 and 8.0 /iin . The extinction law was found 
to be very similar in both fields, and in closer agreement with the obser vations of L96 and 



there fore with models of denser clouds with larger grains. The studies of 



Jiang et al. 



(12003 



20061 ) used data from the ISOGAL survey to estimate the extinction at 7 /xm relative to 
K, and their data showed a good agreement with 105 in a relatively large number of fields 
across the Galactic plane, perhaps suggesting a universal MIR extinction law for clouds 
across the Galactic plane. 

However, one very important aspect to consider is that the regions of the Galactic 
plane studied to date have only moderate extinctions. In particular, the fields studied by 
105 scarcely contained material with very large column densities: one of their fields points 
towards I = 42 deg, which despite coinciding to some extent with with clouds in the Galactic 
Ring, W50 and the Sagittarius Arm, is mostly high above the plane and thus covers 
probably more diffuse material. The other line of sight of 105 crosses the Carina Arm, and 
coincides with a region dense enough to have ongoing star formation, but given the limited 
sensitivity of the 2MASS catalog, they could not obtain reliable colors beyond H — K =1.2, 
and thus their study is limited to zones with extinctions below Ak <2 mag. It would be 
important now to determine if the extinction law remains universal toward regions at much 
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higher extinction levels (2 < Ak < 10), like those that characterize molecular cloud cores 
that form stars. At such regions, extinction can easily exceed the depth probed by 105, and 
the properties of dust might vary significantly because of growth and coagulation of grains. 

In order to test the universality of the extinction law at extreme opacity, we combined 
deep, high resolution NIR (J, H and Kg) observations from ESO-SOFI/NTT and 
ISAAC/ VLT with recently available Spitzer Space Telescope observations in the four IRAC 
MIR bands ([3.6], [4.5], [5.8] and [8.0] /im) of the dark core Barnard 59 (B59, LDN1746), 



which is the most opaque region of the Pipe Nebu 



currently forming a star cluster (IBrooke et al 



a (£=2.9, 6=7.12, d =13O pc) and is 



20061 ). B59 also has the advantage of being 



projected against the highly uniform background of the Galactic bulge, which guarantees 
that any variations of the intrinsic colors of background sources will be insignificant. 
Our data allowed us to measure color excesses for numerous sample of deeply extincted 
stars behind the core otherwise hidden to NIR observations alone. By using the color 
excess of such sources we are able to extend the reddening law between 1.2 to 8.0 /mi to 
unprecedented regimes of extinction. 



2. Observations 



We selected ten SOFI/NTT and two ISAAC/ VLT fields from our large scale NIR 
survey of the Pipe Nebula, which is part of a multi-wavelength project to study the 
primordial conditions of star formation in this nearby, very young molecular cloud. The 
selected fields cover the regions of highest exti nction for B59 and its surroundings, as 



Lombardi et al 



(120061 ). In Figured] we plot 



determined from the 2MASS extinction map of 
the location and extent of the NTT and VLT observations along with the coverage of the 
IRAC mosaic to show the area of coincidence of the surveys. All SOFI/NTT fields were 



observed in H and Kg, and four of them were also observed in J. The ISAAC/VLT fields 



- 7- 



were only observed in H and K$- 



The observations at NTT were carried out during the month of June in 2002 and 
2003, while the VLT observations were done in July, 2003. Raw data were processed using 



standard IRAF pipelines, which were modified from t 
survey of Giant Molecular Clouds. One pipeline (see 



hose used in the NOA Q /FLAMINGOS 



Roman-Zuniga 



20071 ) processes all raw 



frames by subtracting darks and dividing by flat fields, optimizing signal to noise ratios by 

rames. The product 



Levine 



20061), which 



means of a two pass sky subtraction method, and combining reduced 
images from the first pipeline were analyzed by a second pipeline (see 
extracts all possible sources from a given field, applies PSF photometry, calibrates to a zero 
point and calculates an accurate astrometric solution. Both photometry and astrometry 
were calibrated with respect to 2MASS. 

The resultant photometry for the SOFI/NTT fields is reliable {crj t H,K s <0.1 mag) down 
to J=20.4, H=19.2 and #=18.6 mag. The ISAAC/VLT observations, aimed to increase as 
much as possible the detection of sources near the core center, yield reliable photometry 
down to #=22.2 and #5=20.9. 

In the case of the Spitzer observations, we made use of the publicly available pipeline 
catalog of B59 wh ich is part of the Cores to Disks Legacy Project (Data Release 3, S11.0) 



(jEvans et al 
[8.0]. 



20051 ) and includes photometry in the four IRAC bands, [3.6], [4.5], [5.8] and 



3. Analysis 



The relative extinction curve, A x /A Ks (A=[3.6] to [8.0]) for B59 was constructed 
following a similar technique to the one used by 105, except that we calculated the 
relevant color ratios with respect to H — Ks instead of J — K$- We assumed a fixed ratio 
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Ah/Ak = 1.55, which corres ponds to the one d e termi ned by 105 and agrees well with 



the value of 1.56 obtained by 



Rieke fc Lebofsky 



(Il985l ). This choice serves as a simple 



normalization of our law to the one of 105, and any variations in our determinations will be 



relative to their work. With this value, the relative extinction is calculated as: 



A x _ , A H _ x E x . Ks + L 



A Ks A Ks E H _ Ks 

The quotient E\_k s I Eh-k s is simply the slope of the distribution of sources in the 
color-color diagram (A — Ks vs. H — Kg). The slopes were calculated with a weighted 
least square fit applied to restricted color samples. A first restriction limited the allowed 
photometric uncertainties to 0.1 mag in each NIR and MIR band. A second restriction 
assured that only sources with no intrinsic infrared excesses were considered, by separating 
out from the combined catalog every star that had a SED indicative of intrinsic excess (i.e. 



YSOs and galaxies), and keeping only those so urces that fa 



object type categories in the IRAC catalog (see 



Evans et al 



1 in 



2005 



;he 'star' or 'two band' 
table 8). 



As a result of these restrictions, the total number of sources used to calculate the linear 
fits in the J - K vs. H - K s , [3.6] - K s vs. H - K s , [4.5] - K s vs. H - K s , [5.8] - K s 
vs. H - K s and [8.0] - K s vs. H - K s diagrams are 672, 2779, 2776, 1701 and 1170 
respectively. From now on, we will use the nomenclature R[J], R[3.6], #[4.5], #[5.8] and 
i?[8.0] when referring to the color-color diagrams constructed with these samples. 



The least square fits were calculated with a technique similar to that of lLombardi et al. 



( 20061 ): First we divided the H — Ks data set into bins of equal size. The bins have widths 



of 0.056 in R[J], 0.11 mag in #[3.6] and #[4.5], and 0.14 mag in #[5.8] and #[8.0]. These 
bin widths are equivalent to dividing the H — K s data set in bins of constant 'zero excess' 



- 9- 

extinctiorl^l: 0.1 mag for R[J], 0.2 mag for i?[3.6] and _R[4.5], and 0.25 mag for i?[5.8] and 
i?[8.0]. These bin sizes were set by the requirement that a minimum number of sources 
are included in each bin to assure more robust statistics. Then we calculated the median 
and the standard deviation of the two corresponding colors in each bin and applied a least 
squares linear fit with chi-square minimization to the fiducial defined by these median 
points; lastly, we iterated the process by recalculating the fit with only those objects that 
are within 5a of the median points, until the value of the slope converges within 10 -4 . 

In Figures [2] and [3] we show the linear fits applied to samples in the diagrams R[J] 
to i?[8.0]. As mentioned above, the slope of each distribution is the corresponding color 
excess ratio. However, as an independent check we also calcula t ed th e color excess ratios 



Kenyon et al 



(119981) (hereafter KBL), in 



by applying a method similar to the one used by 
which a set of color excesses for each source was calculated by individually subtracting from 
it the colors of each star in a control field. 

Our control field is located approximately two degrees West of Barnard 59 and at the 
same galactic latitude. We also have NIR ESO observations as well as C2D catalog for 
this field. The measured photometric dispersion for sources in these control areas is less 
than 0.08 mag in all the A — K$ colors, which is comparable to the median photometric 
uncertainty of the colors, showing that the bulge field in the background of the Pipe Nebula 
has such a small intrinsic color dispersion that our photometry cannot even resolve it. 
The resultant color-excess data sets are plotted against each other and again, the slopes 
estimated from the least square fits gave directly the ratios of color excess. Comparison 
of the results from the two fitting methods was particularly useful for the R[J] diagram, 
because the number of J detections was small compared to those in the other bands. In the 



2 The 'zero excess' extinction is defined as Ak 3 ,o = 
intrinsic value (H — K)q = 0.0. 



1.78 * (H — Ks), which assumes an 
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case of the color-color method some bins had very few sources and therefore poor statistics, 
which resulted in a least square fit which was also poor, with a slope uncertainty of 0.18. 
When the color-excess method was applied, the statistics per bin became more robust and 
the uncertainty in the slope was reduced to only 0.05. The slopes calculated for the source 
distributions on the color-color and color-excess diagrams, agree within the error in all 
cases, R[J] to #[8.0]. 



4. Results and Discussion 

In Table [1] we list the slopes of the linear fits from the color-color and color-excess 
diagrams for R[J] to i?[8.0]. For comparison, we also list the corresponding slopes for 
the average data points of 105, inferred from the factors A\/Ak quoted on their paper. 
The table also presents the relative extinction ratios Ax/Ak s , which define the extinction 
law determined in this study. The uncertainties listed are those calculated from the 
uncertainties in the slopes of our color-color and color-excess diagrams, and do not include 
the uncertainty of 0.08 in the ratio Ah/Ak s = 1.55 of 105. 

Our distributions, just as those of 105, are well fit with a single slope: we do not see 
any clear evidence for a variation in the slope of the source distributions with increasing 
extinction, even though the background sources we observe through Barnard 59 are much 
more obscured (Ay < 59, A K < 7 mag) than those studied by 105 and other authors (e.g. 
Av <n, A K <2in 105). 

Our study confirms that the extinction law in dense molecular clouds is approximately 
'flat' between 3 and 8 /im. In Figure 3] we show how our ext i nction law compares to the 



'case B' synthetic extinction curves of IWeingartner &: Draind ( 12001 



hereafter WD01) for 



Rv=3.1 (diffuse cloud) and Ry=5.5 (dense cloud). The case B model of WD01 keeps a 
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constant abundance of carbon and silicates, resulting in a distribution of grain sizes that 
favors fewer small grains and a more si gnificant contribution of large grains in dense cloud 



dust, as suggested by observations (e.g. 



Kim &: Martin 



1996). The Ry=5.5 curve is in much 



closer agreement with our data and those of 105 and L96 than the iiV = 3.1 curve. In general 
terms, our data confirm that towards dense molecular cloud material, the MIR extinction 
law departs from the simple extrapolation of the NIR power-law to longer wavelengths. 
The observed opacity in dense clouds appears to become higher and more gray in the MIR, 
with a soft inflexion towards the silicate peak at 10 /im. 

An extinction law that departs from the NIR power-law in such manner could be 
showing how absorption by dust might become more gray in a very dense cloud due 
to anticipated grain g rowth at high column density regimes. In current models (e.g. 



Li fc Greenberg 



20031 ). as a molecular cloud becomes more dense (e.g. when evolving 
towards protostellar collapse), grains are expected to grow via coagulation of smaller grains 
and via the accretion of an organic refractory mantle composed of ices of species depleted 
from the gas phase. In t he first case, there is a net removal of small grains (0.1 /xm) and 



PAHs in the dense cloud ( Chokshi et al. 



19931 ). which results in a significant contribution to 



extinction by particles with larger cross sections. In the second case, the contribution of an 
accreted ice mantle to the total extinction at any given wavelength is expected to be less 
i mportant be cause the mantles may not increase the grain sizes by more than about 18% 



( IDraine 



19851 ). In any event, although our observations may suggest that the dust grains in 



a dense cloud are larger than those expected for diffuse gas, we find no significant evidence 
for additional grain growth as a function of increasing extinction, and presumably density, 
within B59 itself. 

The extinction law we determined towards B59 appears to be different from the one 
inferred from the SWS observations of the Orion Molecular Cloud by 



Bertoldi et al 



(1999) 
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clouds (e.g. 



Li fc Draine 



and iRosenthal et al.l (120000 . w hich is similar to the one expected from models of diffuse 



20011 ). However it is not clear that the OMC is a good example of 
a diffuse cloud: even though the Orion cloud is located at a high galactic latitude, meaning 
a very small contamination by foreground dense clouds, one has to consider that much of 
the dusty material observed toward the OM C probably arises in t he HII region and in dense 



clouds behind it. Indeed, as pointed out by lBertoldi et al. 



(119991 ) their "effective extinction 



may not necessarily follow an average inters tellar extinction l aw, si nce the emitting and 



absorbing gas may be mixed". In addition, 



Rosenthal et al. 



(120001 ) warned that there 



are large uncertainties in their calculations, arising mostly from their calibration and the 
presence of foreground H2 fluorescence, but that the minimum of opacity at 6.5 /im (which 
implies the continuation of the NIR power law) they observe could still be real. 

The slopes of the distributions of B59 sources in our color-color diagrams are less steep 
than those inferred from the relative extinction ratios A\/Ak s calculated by 105 towards 
two distinct fields of the Galactic plane (see figures [2] and [3]) . As a result of this, the 
extinction law, A\/ Ak s as a function of A, towards B59 is also less steep than the one 



determined by 105. 



Chiar fc Thielend (120061 ) derived an extinction curve in the form of a 



polynomial fit to the combined data points of 105 and L96, and we compare this curve to 
our data points in Figure HI The difference we observe is modest, especially considering the 
large uncertainties in the relative extinction ratios of 105 and the data points of L96, but it 
could indicate that the extinction law in the NIR and MIR may vary slightly as a function 
of environment. In particular, 105 and L96 observations not only sample different lines of 
sight but also a not well constrained mixture of diffuse and dense cloud material, while our 
observations only sample dense cloud dust. 

Our study was successful at combining NIR and MIR photometry to derive the 
extinction law through a highly opaque core in a molecular cloud. Barnard 59 is a good 
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example of a cloud so obscured that NIR observations alone cannot penetrate it, and 
thus to detect and determine the colors of stars in its background, MIR observations are 
required. This situation is common to the centers of other very dense pre-stellar and young 
star forming cores, and even relatively thin IRDCs. In our experiment, the extinction 
law is measured specifically at opacity levels similar to those of such regions, and can 
be used in combination with deep NIR and MIR observations to probe highly obscured 
clouds with more confidence than if one had to use an extinction law extrapolated from 
measurements at lower opacities. It is thus possible to construct very deep extinction 



maps u sing Spitzer data, as we will show in an upcoming paper (IRoman-Zuniga et al. 



in prepl ). where we use Spitzer observations and our derived extinction law to construct an 
extinction map for B59 that reaches a depth close to Ax s = 10 mag, and resolves the most 
opaque regions near its center. Detailed knowledge o f the inner structure of such a core can 



provide critica 



Teixeira et al 



2004 



information about its physical state (lAlves et al 



Harvey et al 



2001 



2005 ), and its stability against gravitational collapse (IShu 



2003|). 



Kandori et al 



1977 



200 



Lada et a] 
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Fig. 1. — Common observations of B59. The large trapezoidal area denotes the coverage 
of IRAC. The irregularly shaped medium size area delineates the span of the NTT obser- 
vations and the small region at the center indicates the fields observed with the VLT. The 
common, shaded region is the one used in this study. The halftone and contours in the 
ba ckground indicat e extin ction Ay from 5.0 to 25.0 mag, in steps of 2.5 mag, as determined 



by 



Lombardi et al. 



feoOJ). 
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Fig. 2. — J — K vs. H — Ks diagrams for sources with counterparts in both the Spitzer and 
ESO catalogs for the region of Barnard 59. Top: color-color diagram, Bottom: color-excess 
diagram (control sources are from an extinction-free field at the same galactic latitude as 
B59). In both plots gray points are ESO sources, the dotted line is the fit to the median 
points indicated by the large solid dots and the dashed line represents the median slope for 
sources in the fields studied by 105 -which has been shifted arbitrarily in the vertical direction 
to adjust to the ordinate of our fit. 
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Fig. 3. — Color-color diagrams for combined ESO/IRAC colors. Top left: [3.6] — K vs. 
H - K s (R2); top right: [4.5] - K vs. H — K s (R3); bottom left: [5.8] — K vs. H - K s 
(R4); bottom right: [3.6] — K vs. H — Ks (R5). The layout of symbols and lines is the same 
as in Fig. El 
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Fig. 4. — The extinction in Barnard 59 as determined with our analysis. Open circles 
represent the relative extinction ratios A\/Ak s calculated directly from the fits to points 
in color-color diagrams. Solid circles symbols are extinction ratios calculated from the fits 
to the distributi ons of color excesses (KB L method). The polynomial fit to data points of 



105 and L96 by 



Chiar fc Thielend (120061 ) (CT06) is shown with a solid line. The 'case B' 



synthetic extinction curves of WD01 for Ry=3A and 5.5 are plotted with dot-dashed and 
dashed lines respectively. 
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Table 1. The Extinction Law between 1.2 and 8.0 /jm in Barnard 59 



A a E(X-K S )/E(H-K S ) A x /A Ks 

[fim] 105 b C-C c C-E d C-C C-E 



1.240 2.73 

1.664 1.00 

2.164 0.00 

3.545 -0.80 

4.442 -1.04 

5.675 -1.04 

7.760 -1.04 



2.36±0.18 
1.00±0.00 
O.OOiO.OO 
-0.66±0.02 
-0.83±0.03 
-0.97±0.02 
-0.96±0.02 



2.52±0.05 
1.00±0.00 
O.OOiO.OO 
-0.69±0.03 
-0.89±0.04 
-0.98±0.05 
-0.97±0.04 



2.299±0.530 
1.550±0.080 
1.000±0.000 
0.618±0.077 
0.525±0.063 
0.462±0.055 
0.455±0.059 



2.389±0.045 
1.550±0.080 
LOOOiO.OOO 
0.619±0.029 
0.512±0.022 
0.459±0.024 
0.465±0.017 



a Isophotal wavelengths of 2MASS and IRAC convolved with a K2 III 
star (see 105). 

b Slopes of 105 inferred from the A\j A^ s ratios of their Table 1. 
c Slopes of the distributions of source colors. 
d Slopes of the distributions of source color-excesses. 



